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Abstract 

We compute the Coulomb effects produced by an expanding, highly charged 
fireball on the momentum distribution of pions. We compare our results to data 
on Au+Au at 11.6 A GeV from E866 at the BNL AGS and to data on Pb+Pb at 
158 A GeV from NA44 at the CERN SPS. We conclude that the distortion of the 
spectra at low transverse momentum and mid-rapidity can be explained in both 
experiments by the effect of the large amount of participating charge in the central 
rapidity region. By adjusting the fireball expansion velocity to match the average 
transverse momentum of protons, we find a best fit when the fireball radius is about 
10 fm, as determined by the moment when the pions undergo their last scattering. 
This value is common to both the AGS and CERN experiments. 
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The purpose of colliding heavy nuclei at relativistic and ultrarelativistic energies is to 
produce high density, high temperature matter, either hadronic or quark-gluon plasma. 
One way of inferring the properties of matter under these extreme conditions is to compare 
the single particle momentum distributions of secondaries to those obtained in collisions 
of lighter systems at similar beam energies. An example is provided by the qualitative 
differences found between the transverse momentum distributions from heavy-ion reac- 
tions and those from p-p collisions. In the former, the bulk features of the spectrum can 
successfully be explained by the large amount of secondary scattering, potentially leading 
to the onset of hydrodynamical behavior jj]]. 

Another feature to account for in the collision of heavy systems is the presence of 
a large amount of electric charge that will influence the dynamics of secondary charged 
particles. Due to the long-range nature of the Coulomb field, the spectrum of charged 
particles can still be distorted even after freeze-out. For central collisions this Coulomb 
effect can be more significant in the case of strong stopping when the participant charge 
in the central rapidity region is an important fraction of the initial charge, particularly 
on the low momentum particles. 

In a recent paper || two of us explored the influence of the time dependent electric field 
produced by an expanding, highly charged fireball on the low energy part of the spectrum 
of charged test particles, taken to be kaons [[| . The analysis was based on the solution to 
Vlasov's equation assuming that the test particles' dynamics was nonrelativistic. In this 
paper we extend that analysis, treating the test particle relativistically. We also apply 
the analysis to the transverse momentum distribution of charged pions and compare the 
results to recent measurements of Au+Au collisions by the E866 collaboration at the BNL 
AGS and to recent measurements of Pb+Pb collisions by the NA44 collaboration at the 
CERN SPS. Relativistic calculations for pions and kaons have also been done in somewhat 
different models by Barz et al. Q. 

We briefly recall the basics of the model used in Ref. [pj and how it is modified for 
relativity. A uniformly charged sphere which has a total charge Ze and whose radius R 
increases linearly with time t from a value Rq at time to at a constant surface speed v s 
produces an electric potential 



In the center-of-mass frame of the fireball the charge moves radially outwards, hence 
there is no preferred direction and consequently the magnetic field produced by this 
moving charge configuration vanishes. The fireball parameters are related by Rq = v s to- 
If / ± (r,p,t) represents the ±e test particle phase space distribution then, when ignoring 
particle collisions after decoupling, its dynamics is governed by Vlasov's equation. 




r > R = v s t 
r < R = v s t 



(1) 



+ |_. Vr ±eE(r,t)-V p / ± (r,p,t) = 0, 



(2) 
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where E p = y/p 2 + m 2 , m is the meson's mass and 



eE(r, t) = -eV r V(r, t) = { <f>lf} J ' r ^ = ^ (3) 
v ' y ' I (t s vf/4r 6 ) r , r > R = v s t y ' 

is the time-dependent electric field corresponding to the potential V(r,t). Here we have 
defined the characteristic time t s by 

ts = 3 , 4 

and we work in units where c = 1. The set of Eqs. d3) — C@) are relativistically correct. 
Notice that there is no problem with retardation effects since what is specified in these 
equations is the electric field and not the charge density. Of course, the relationship 
between the charge density and the eletric field gets modified due to the finite speed of 
light, but we are not interested in that relationship here. The only difference between the 
equations to be solved here and those solved in Ref. involve the energy E p in Eq. (|2|) 
which was approximated by m + p 2 /2m in the nonrelativistic case. 

The solution to Eq. (0) is found by the method of characteristics. This involves solving 
the classical equations of motion and using the solutions to evolve the initial distribution, 
assumed here to be thermal 

/ ± (r, p, t ) = exp {- (E p ± V{r, t )) /T} (5) 

forward in time. The equations of motion are: 

dr p 

dt ~ ~W P 

f = ±eE(r,t). (6) 

The pion's asymptotic momentum is calculated numerically by a 6'th order Runge-Kutta 
method with adaptive step sizes from a set of initial phase-space positions. The final 
momentum distribution is the result of computing the trajectories for many initial phase 
space points. The initial radial position was incremented in N r =50 steps with spacing 
At/Rq = 0.02. The initial momentum was incremented in N p =600 steps with spacing Ap 
= 1 MeV. The cosine of the angle between the initial position and momentum vectors 
was incremented in N 2 =100 steps of size A cos (9 = 0.02. Hence the total number of 
trajectories computed was 3 x 10 6 for each set of initial conditions. 

For central Pb+Pb collisions at 158 A GeV at the SPS, for which NA44 and NA49 
have reported a significant amount of stopping |5|, [|, the number of protons per unit 
rapidity in the central region is on the order of 30. If we consider that the rapidity region 
spanned by the fireball is between 1 and 5 units, then we take as the effective fireball's 
charge Z = 120. 
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To estimate the surface velocity v s we note that the average transverse velocities of 
pions and protons in Pb+Pb collisions at the SPS are greater than 0.6 0. Therefore, to 
represent the electric field it is reasonable to take the maximum velocity of charges to be 
greater than that; we shall use v s = 0.8 and v s = 1. 




Figure 1: The transverse mass distributions at mid-rapidity for tt + and 7r . The solid 



lines are exponentials with an inverse slope of 110 MeV. The data are from NA44 [14 



The model we are using is spherically symmetric. It has long been known @ that the 
momentum distributions at the SPS are somewhat forward-backward peaked, even for 
central Pb+Pb collisions ||. A spherically symmetric model is, of course, not essential 
to the basic physics and can be relaxed at the expense of additional computing time. 
Nevertheless, we can provide the following arguments in favor of its use. First, we will be 
comparing with the transverse momentum distributions at mid-rapidity where the impact 
of spherical asymmetry should be less important than near the fragmentation regions. Sec- 
ond, pion interferometry of central Au+Au collisions at the AGS || and of central Pb+Pb 
collisions at the SPS [|10j both yield comparable values for the transverse and longitudinal 
radii at the time of pion freezeout or strong decoupling. These are about twice the radius 
of a cold gold or lead nucleus. Third, we shall show later that the transverse surface of 
the fireball must expand outwards with a speed more than 90% that of light in order to 
reproduce the average proton transverse momentum. Since the longitudinal surface of the 
fireball cannot travel faster than the speed of light, this means that in velocity space the 
fireball is nearly symmetric. These phenomena, although not yet measured at that time, 
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were already known to Landau ||TTJ. His model envisioned the two nuclei stopping each 
other within one Lorentz contracted nuclear diameter, then undergoing a one dimensional 
expansion along the beam axis. After a time of order (nuclear radius) / c the longitudinal 
size of the fireball had become comparable to the transverse size and, furthermore, the 
rarefaction wave would have caused transverse expansion to set in. As a rough approx- 
imation one therefore could imagine the original one dimensional expansion to become 
three dimensional after that time. This model was later developed by others, including 



Cooper, Frye and Schonberg |12[ . The essential insight from Landau's model is not the 
degree of stopping but rather the point that significant transverse expansion sets in after 
the longitudinal and transverse radii become comparable in magnitude. Thereafter, from 
the point of view of a distant observer, the expansion is not as asymmetric as one might 
originally think. 




Figure 2: The ratio tt + /tt~ vs. ttit — m at the SPS. The data is from ||. Top panel: The 
curve corresponds to the Coulomb calculation described in the text with surface velocity 
0.8 and initial radius 13 fm. Bottom panel: Same as top panel but with the parameters 
1 and 8 fm, respectively. 

The collaboration NA44 has reported a suppression of the tt + /tt~ ratio as — > in 
Pb+Pb collisions, presumably due to the Coulomb field of the fireball. Only a very small, 
or no, suppression is observed in S+Pb collisions. To model the primordial distribution 
we use an exponential parametrization of dN/p 2 dp. Fig. 1 shows this representation [|L| 



in comparison to the transverse mass distributions of positive and negative pions [14 with 
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a value of T e g = 110 MeV. This common fit (apart from the absolute normalization) is 
then used to model the initial pion distribution, and the Coulomb computation we do is a 
perturbative correction to it. In essence we sprinkle the pions around at time to according 
to a relativistic Boltzmann distribution with an effective temperature of 110 MeV and 
solve the Vlasov equation for subsequent time. This assumes that all resonances that 
can decay into charged pions have already done so, hence their effects are contained in 
this parametrization. Even most of the relatively long-lived omega-mesons should have 
decayed by the this late freeezeout time and relatively low temperature. 

The procedure we follow is to fix v s at either 0.8 or 1 and then vary Rq to get the 
best fit to the ratio tt~ /tt + as a function of my — m in the range from to 500 MeV. 
The data's bin size is 10 MeV. The data are shown in Fig. 2 and compared to the 
results of calculations with v s = 0.8, Ro = 13 fm and v s = 1, Ro = 8 fm. There is no 
visual difference between the calculations with the two parameter sets. The reason the 
final Coulomb corrected ratio is nearly identical for these two sets of parameters is easily 
understood. According to the study in Ref. 0, the longer the test charge stays in the 
electric field of the fireball the more it is modified from its initial form. Thus, for the 
more slowly expanding fireball (v s = 0.8) the radius must be larger (13 vs. 8 fm) in order 
to reduce the magnitude of the electric field so as to represent the data. 
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Figure 3: Chi-squared vs. radius for fixed values of the surface velocity for pions for 
central Pb+Pb collisions at the SPS. The curves represent polynomials which pass through 
the values of Ro for which computations were done. 
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In Fig. 3 we plot the chi-squared for the 7T + /tt~ ratio as a function of radius Rq for the 
surface velocity v s fixed at 0.8 and 1. The minimum is achieved at about 13 and 8.3 fm, 
respectively. The difference in the chi-squared for the two values of v s is not significant 
enough to make a judgement on which gives the better fireball physics. The reason for 
the close coupling of the two parameters, v s and Rq, is explained by recognizing that it 
is the ratio of the characteristic expansion time to to the characteristic Coulomb time t s 
which is the relevant quantity. Since t /t s oc v 2 s Rq we would expect comparable fits to the 
data will be obtained when Rq varies inversely with the square of v s ] this is borne out 
by our explicit numerical calculations. We remark that a different parametrization of the 
primordial pion distribution, such as a two temperature fit, might well lead to an even 
better representation of the data. 
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Figure 4: Chi-squared vs. radius for a surface velocity of 0.914 for pions for central 
Au+Au collisions at the AGS. The curves represent polynomials which pass through the 
values of Rq for which computations were done. 

The immediate question now is how to determine the appropriate value of v s . We 
shall do this by matching the average transverse momentum of the protons as inferred 
from our assumption of a uniformly expanding sphere with the measured one. The model 
gives 

, . / d 3 r p(r,t) smOmpv/^/l — v 2 



where p(r,t) = po(t /t) 3 Q(v s t — r) and v = r/t. The result is: 



2- (2 + v 2 s )Jl-vl 
^ = n 7^3 m P- ( 8 ) 

a 

In the limits v s —* and f s = 1 one gets (py) = 37rmpt> s /16 and (pr) = nmp/2, 
respectively. Note that v s should not be interpreted as a hydrodynamic flow velocity. 
Rather, it embodies the combined effects of hydrodynamic flow and thermal motion of 
the net charge carriers, mainly protons. Both NA44 ]TJj] and NA49 |T5| have reported the 
transverse mass distribution in central Pb+Pb collisions at midrapidity to be dN/pxalpx oc 
exp(— ttit/Tp) with Tp = 290 MeV. This corresponds to an average transverse momentum 
of 825 MeV. The value v s = 0.916 gives the best matching of the model to the proton 
spectra. This finally pins down the freezeout radius for pions to be 10 fm. 
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Figure 5: The ratio 7r /tt + vs. 
radius at pion freezeout is 10 fm. 



m at the AGS. The data is from 16 . The fireball 



We now repeat the above analysis for the BNL experiment E866 [16|. This experiment 
measured protons and positive and negative pions for Au+Au collisions at 11.6 A GeV. 
We focus on the most central events, those with the highest multiplicity and transverse 
energy. The average proton transverse momentum is 820 MeV yielding the surface speed 
v s = 0.914. The number of participating protons maybe be estimated from the data to be 
Z = 116. The slopes of the ir + and ir~ are the same to within several MeV; their average 
is 111 MeV. All of these quantities are remarkably similar to those in central Pb collisions 
at the much higher energy of the SPS. 
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The experiment reported the ratio n~ /n + for a central window of rapidity. We cal- 
culated the ratio with two parameters: the radius at pion freezeout Rq and the absolute 
normalization of the individual spectra. The chi-square per degree of freedom is shown in 
Fig. 4 as a function of Rq. It has a minimum of 0.53 near 10 fm, the same as Pb collisions 
at the much higher SPS energy! The lower value of chi-square in this case as compared 
to the SPS data may be due to the better approximation of spherical symmetry. The 
calculated ratio is compared to the data in Fig. 5 and its appearance is quite satisfactory. 

In conclusion, we have demonstrated that the suppression of the ratio tt + /tt~ in central 
Pb+Pb collisions at the SPS as observed by NA44 and in central Au+Au collisions at 
the AGS as observed by E866 can be quantitatively understood as a Coulomb effect. The 
large electric fields generated by the expanding fireballs seem to distort the pion ratios 
at small transverse momenta and midrapidity in accordance with these experiments. In 
principle these ratios should provide a very good measure of the size of the fireball at 
the time pions cease collisions and begin to free-stream. From experiments |TJ one may 
estimate that about 2500 baryons and mesons are produced with rapidities within 2 units 
of mid-rapidity in the central Pb+Pb collisions. Associating these with a fireball of radius 
10 fm means an average density in the fireball reference frame of 0.60/fm 3 ; this is 4 times 
nuclear density! Taking into account relativistic length contraction the local density is 
smaller by a factor 



For v s = 0.916 this implies a local density of 0.42/fm 3 , still quite considerable. This 
implies that at the moment that free-streaming of mesons begins the hadronic density is 
still interestingly high. 

To make further progress it is clear that more accurate data are called for but, more 
importantly, one should also compute the influence of the Coulomb force within an event 
simulator which takes into account a more accurate description of the space-time evolution 
of the collision. 
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